Search results

The detailed flow behaviour around a four—hole Cobra Pitot pressure probe, developed by the Commonwealth Scientific and Industrial Research Organization, Australia, (CSIRO), to determine the pressure and the velocity components in three dimensional single‐phase/multi‐phase fluid flow, is investigated computationally. The incompressible steady state Navier—Stokes equations are solved numerically using a general purpose computational fluid dynamics (CFD) code developed at CANCES. Computational results are presented for representative probe pitch and yaw angles at a Reynolds number = 2 × 10³, emphasising the pressure distribution and flow separation patterns on the probe tip adjacent to the pressure ports. Quantitative comparison of the computational simulation to experimental results is done by comparing experimental calibration data to numerically computed pressure responses. The topological features of the near tip flow behaviour are visualised using critical point concepts and three dimensional streamlines. Additional qualitative comparison to experiment is discussed using data from a preliminary experimental investigation using surface oil film visualisation techniques, where available. Conclusions are drawn concerning the near tip flow behaviour, the good level of agreement between the numerical results and experimental data and the effectiveness of using a computational analysis to provide accurate detail useful for engineering design purposes.

An alternative algorithm has been developed for computing the behaviour of flows within arbitrary ducts and channels. This technique requires a small number of downstream marches in the primary flow direction, employing, on each march, numerically efficient procedures originally developed for a single sweep non‐elliptic flow solver. The multiple sweeps allow the capture of effects such as upstream pressure influences and streamwise recirculation. The energy equation is also solved to allow for varying heat transfer between the fluid and the boundary walls. The numerical work is further complicated by considering flows within turning sections of ducts which demonstrate large transverse velocities and consequent distortion of the primary flow. The computations are validated by comparison with a number of fluid/heat transfer experiments. The majority of these are taken from studies of turning flows within circular arc ducts which display the various pressure and transverse flow phenomena for which this new algorithm was initially developed to represent.

Numerical predictions of laminar and turbulent fluid flow and heat transfer around staggered and in‐line tube banks are shown to agree closely with seven experimental test cases. The steady state Reynolds‐averaged Navier‐Stokes equations are discretised by means of a cell‐centred finite‐volume algorithm. Two‐dimensional results include velocity vectors and streamlines, surface shear stresses, pressure coefficient distributions, temperature contours, local Nusselt number distributions and average convective heat transfer coefficients, and indicate very good agreement with experimental data. It is found that a relatively fine grid is required to be able to predict the surface heat transfer behaviour accurately. Also, three‐dimensional simulations are shown, which are physically consistent. The numerical procedure presented here is robust, accurate and time efficient, making it suitable as a design tool for tube banks in heat exchangers.

A modelling approach of gas solid flow, considering different physical phenomenon such as fluid turbulence, particle turbulence and interparticle collision effects are presented. The approach is based on the two‐fluid model formulation where both phases are treated as continuum. This implies that the gas phase as well as the particle phase are weighted by their separate volumetric fractions. According to the experimental results and numerical simulations, the inter‐particle collision possesses a significant influence of turbulence level on particle transport properties in gas solid turbulent flow even for dispersed phase volume fraction (α<0.01). Comparisons in predictions have been depicted with inclusion of interparticle collision effect in the equation of particle turbulent kinetic energy and with exclusion of this effect. Experimental research has been conducted in a thermal power plant depicting higher erosion resistance of noncircular square sectioned coal pipe bends in comparison with those with circular cross section, the salient features of the experimental work are presented in this paper. Experiments have been conducted to determine, pressure drop in straight and curved portions of conduits conveying air coal mixtures in a thermal power plant. Validation of this experimental data with numerical predictions have been presented.

The Least Restrictive Environment (LRE) requirements of the Individuals with Disabilities Education Act address factors to consider in educating students with and without disabilities together to the maximum extent appropriate. This chapter is designed to examine the origins and evolving interpretations of the LRE concept in special education policy and practice. Discussion traces the evolution of the concept as a legal principle, and reviews its application to educational strategies for students with learning and behavioral disabilities in contemporary schools. In conclusion, the future of the LRE concept is addressed in light of competing policies promoting presumptive inclusive education, and publicly funded school choice programs promoting greater involvement of parents in choosing where their children with and without disabilities should be educated.

Detailed results of numerical calculations of transient, 2D incompressible flow around and in the wake of a square prism at Re = 100, 200 and 500 are presented. An implicit finite‐difference operator‐splitting method, a version of the known SIMPLEC‐like method on a staggered grid, is described. Appropriate theoretical results are presented. The method has second‐order accuracy in space, conserving mass, momentum and kinetic energy. A new modification of the multigrid method is employed to solve the elliptic pressure problem. Calculations are performed on a sequence of spatial grids with up to 401 × 321 grid points, at sequentially halved time steps to ensure grid‐independent results. Three types of flow are shown to exist at Re = 500: a steady‐state unstable flow and two which are transient, fully periodic and asymmetric about the centre line but mirror symmetric to each other. Discrete frequency spectra of drag and lift coefficients are presented.

Variable density flows play an important role in many technological devices and natural phenomena. The purpose of this paper is to develop a robust and accurate method for low Mach number flows with large density and temperature variations.

Low Mach number approximation approach is used in the paper combined with a predictor-corrector method and accurate compact scheme of fourth and sixth order. A novel algorithm is formulated for the projection method in which the boundary conditions for the pressure are implemented in such a way that the continuity equation is fulfilled everywhere in the computational domain, including the boundary nodes.

It is shown that proposed implementation of the boundary conditions considerably improves a solution accuracy. Assessment of the accuracy was performed based on the constant density Burggraf flow and for two benchmark cases for the natural convection problems: steady flow in a square cavity and unsteady flow in a tall cavity. In all the cases the results agree very well with exemplary solutions.

A staggered or half-staggered grid arrangement is usually used for the projection method for both constant and low Mach number flows. The staggered approach ensures stability and strong pressure-velocity coupling. In the paper a high-order compact method has been implemented in the framework of low Mach number approximation on collocated meshes. The resulting algorithm is accurate, robust for large density variations and is almost free from the pressure oscillations.

The NIRVANA project is concerned with the development of a nonoscillatory, integrally reconstructed, volume‐averaged numerical advection scheme. The conservative, flux‐based finite‐volume algorithm is built on an explicit, single‐step, forward‐in‐time update of the cell‐average variable, without restrictions on the size of the time‐step. There are similarities with semi‐Lagrangian schemes; a major difference is the introduction of a discrete integral variable, guaranteeing conservation. The crucial step is the interpolation of this variable, which is used in the calculation of the fluxes; the (analytic) derivative of the interpolant then gives sub‐cell behaviour of the advected variable. In this paper, basic principles are described, using the simplest possible conditions: pure one‐dimensional advection at constant velocity on a uniform grid. Piecewise Nth‐degree polynomial interpolation of the discrete integral variable leads to an Nth‐order advection scheme, in both space and time. Nonoscillatory results correspond to convexity preservation in the integrated variable, leading naturally to a large‐Δt generalisation of the universal limited. More restrictive TVD constraints are also extended to large Δt. Automatic compressive enhancement of step‐like profiles can be achieved without exciting “stair‐casing”. One‐dimensional simulations are shown for a number of different interpolations. In particular, convexity‐limited cubic‐spline and higher‐order polynomial schemes give very sharp, nonoscillatory results at any Courant number, without clipping of extrema. Some practical generalisations are briefly discussed.

A finite element based numerical model is employed to obtain isothermal and heat transfer predictions for the case of turbulent flow with a decaying swirl component in a stationary circular pipe. An assessment is made on the quality of predictions based on the choice of turbulence modelling technique adopted to close the governing equations. In the present work the one‐equation, two‐equation and algebraic Reynolds stress turbulence models are employed. For the confined flow problem investigated, accurate prediction of the near‐wall conditions is essential. This is particularly the case for confined swirling flow where the variation of variables near the wall is often somewhat greater than encountered in pure axial flow. A finite element based near‐wall model is employed as an alternative to conventional techniques such as the use of the standard logarithmic functions. Of significance is the fact that flow predictions based on the use of the unidimensional finite element techniques are closer to experiment compared to the wall function based solutions for a given turbulence model. As expected, improvements in the flow predictions directly contribute to improved simulation of the thermal aspects of the problem.

The purpose of this paper is to use the polynomial differential quadrature method (PDQM) to find the numerical solutions of some Burgers'‐type nonlinear partial differential equations.

The PDQM changed the nonlinear partial differential equations into a system of nonlinear ordinary differential equations (ODEs). The obtained system of ODEs is solved by Runge‐Kutta fourth order method.

Numerical results for the nonlinear evolution equations such as 1D Burgers', coupled Burgers', 2D Burgers' and system of 2D Burgers' equations are obtained by applying PDQM. The numerical results are found to be in good agreement with the exact solutions.

A comparison is made with those which are already available in the literature and the present numerical schemes are found give better solutions. The strong point of these schemes is that they are easy to apply, even in two‐dimensional nonlinear problems.